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Abstract 

To examine the applicability of the time-dependent density-functional the- 
ory (TDDFT) for treating the electron-nucleus coupling in excited states, 
we calculate the strength distribution associated with the tt-tt* transition in 
ethylene. The observed optical transition strength at 7-8.5 eV region shows a 
complex structure arising from coupling to C-C stretch motion, to torsional 
motion, and to Rydberg excitations. The mean energy of the observed peak 
is reproduced to about 0.2 eV accuracy by the TDDFT in the local density 
approximation (LDA). The reflection approximation is used to calculate the 
peak broadening. Roughly half of the broadening can be attributed to the 
fluctuation in the C-C coordinate. The asymmetry in the line shape is also 
qualitatively reproduced by the C-C coordinate fluctuation. We find, in agree- 
ment with other theoretical studies, that the torsional motion is responsible 
for the progression of weak transition strength extending from the peak down 
to about 6 eV. The LDA reproduces the strength in this region to about factor 
of 3. We conclude that the TDDFT is rather promising for calculating the 
electron nucleus coupling at short times. 
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I. INTRODUCTION 



The time-dependent density functional theory (TDDFT) offers a compromise between 
economy and accuracy for the calculation of electronic excitations in atoms, molecules and 
condensed matter. The electronic excitation energies are typically reproduced by a few 
tenths of an eV for states with large oscillator strengths, and their strengths are typically 
reproduced within 25% accuracy [1,2]. We are now interested in seeing how much further 
the theory can be applied at this level of accuracy and, in particular, how well it works in 
describing the coupling of electronic excitations to the nuclear degrees of freedom. 

In a previous publication, we examined this question taking the benzene molecule for 
a case study [3]. There we found that the TDDFT worked rather well at describing the 
nominally forbidden transitions, reproducing oscillator strength over 3 orders of magnitude 
with errors at the 30% level. Besides giving transition strength to the symmetry-forbidden 
transitions, the electron-nucleus coupling induces a vibrational splitting of each electronic 
line. While the detailed description requires the Franck-Condon factors of the individual 
vibrational states, the overall widths can be calculated with much simpler theory. For 
benzene, we found that the overall widths associated with the vibrational splittings were 
reproduced within 25% accuracy. 

The ethylene is an even simpler molecule whose absorption spectrum has been subject 
to much theoretical study by the methods of quantum chemistry. The observed photoab- 
sorption spectrum in ethylene shows a strong, structured peak in the region of 7-8.5 eV 
excitation that may be largely attributed to the tt-tt* transition. The studies of ref. [4] and 
[5] present detailed calculations of the vibronic couplings, but they do not make a direct 
comparison to experiment. We also note the calculations by Ben-Nun and Martinez [6] 
in which the subsequent nuclear motion is treated quantum mechanically. The goal of the 
present article is less ambitious: we aim to see how well the density functional theory (DFT) 
works for electron-nucleus coupling at short times, before the system responds dynamically. 
We will focus on the n-ir* transition, which is of general interest for the optical properties 
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of conjugated carbon systems. 

The basic properties of the transition, besides its mean energy, are its width (about 1 
eV FWHM) and its oscillator strength, measured to be / = 0.42 [7]. A feature of particular 
interest is the long progression of strength on the low-energy side of the peak, extending 
down to about 6 eV. We shall also apply the TDDFT to this feature and directly compare 
the theoretical and the experimental strength in the region of 6-7 eV. 



II. CALCULATIONAL METHOD 

We use the DFT to calculate the ground potential energy surface (PES) as well as the 
excited PES. In the present work, we use a simple density functional based on the local 
density approximation (LDA) [8]. 

Denoting the nuclear coordinates by Q, the ground PES E gs (Q) is obtained in the static 
theory by minimizing the energy functional of the orbital variables. The optimized fa satisfy 
the Kohn-Sham equations 

+ ^fa = £i<Pi- (1) 
2m on 

Here £ is the energy functional including Coulomb interactions but not including the elec- 
trons' orbital kinetic energy. 

The excited states are calculated by solving the linearized equations of the TDDFT. 
These equations are very close to the RPA equations, differing only by the presence of an 
exchange-correlation term in the interaction. For each ground state orbital fa, there are two 
excited wave function, <ftfand<f)~ . The equations they satisfy are 

V 2 58 5 2 S 

with the normalization condition, (fai\fat) — {4>i\(t>i) — 1- The physical quantities extracted 
from the solutions are the eigenvalues uj, representing vertical excitation energies, and the 
transition densities 5n = J2i fai^t + which represent the matrix elements of the density 
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operator between the ground and the excited states. The excited PES is given by the sum 
of the ground PES and the excitation energy u, 

E ex (Q)=E gs (Q)+iu(Q). (2) 

To describe the oscillator strength distribution, we treat nuclear motion classically except 
for the zero-point motion on the ground PES. This is called the reflection approximation 
[9,10]. Also, we neglect couplings between different vibrational degrees of freedom and treat 
the fluctuations in various coordinates Q independently. The reflection formula for the 
strength function is 

where the coordinate Q in this expression is obtained by inverting the equation E = E ex (Q). 
In the formula, ^(Q) denotes the vibrational nuclear wave function in the ground PES, 
normalized as / \ ty\ 2 dQ = 1. Also f(Q) denotes the oscillator strength at a fixed Q between 
the ground and the excited state. In ref. [3], we used an even simpler approximation, namely 
taking the vibrations to be harmaonic and treating other factors in eq. (3) as constant. Then 
the width is Gaussian with an energy variance oe given by 

dE P 



ue = 



°Qi (4) 



o 



dQ 

where <tq is the variance of Q in the ground state vibrational wave function. For a Gaussian 
profile, the variance is related to the width T by T = 2.3<7e- 

Our computational method is quite different from usual quantum chemistry methods 
in that we employ a three-dimensional Cartesian grid representation for the orbital wave 
functions. Only valence electrons are explicitly included in the functional. More tightly 
bound electrons are taken into account by pseudopotentials constructed with a commonly 
used prescription [11,12]. Further details on the computation may be found in our previous 
publications, e.g. [1]. The important numerical parameters in the computations are the 
mesh spacing Ax and the shape and size of the volume on which the grid is defined. In our 



previous work on conjugated carbon systems, we used a grid spacing of Ax = 0.3 A. With 
that mesh, orbital energies are converged to better than 0.1 eV. One can use the same coarse 
mesh to compute the coupling to the vibrational degrees of freedom, provided the ground 
PES is independently known. In ref. [3] we only considered small amplitudes of nuclear 
motion, obtaining the necessary information from empirical data on the vibrations. Here we 
want to consider larger amplitudes of motion, beyond the harmonic regime. This requires a 
more accurate calculation of the ground PES, obtainable with a finer mesh. 

III. RESULTS 
A. Ground PES 

We first consider the ground PES. To obtain converged results for equilibrium nuclear 
positions and vibrational frequencies, we found it necessary to use mesh spacing of Ax = 0.2 
A. The orbital wave functions are defined on the grid points that lie within a sphere. We 
found that a sphere of radius R = 5 A is adequate to obtain converged ground orbital ener- 
gies. Some properties related to the ground PES are shown in Table I. The C-C equilibrium 
distance is reproduced to 1% accuracy. The curvature of the potential energy surface for the 
two important modes, the C-C stretch and the torsion, control the respective vibrational 
frequencies. As may be seen from the Table, the empirical vibrational frequencies [13] are 
reproduced to about 5% accuracy. We also quote the results of the ab initio calculations of 
ref. [14] and [5], which show more or less the same accuracy. We conclude that the quality 
of the DFT results are encouraging to proceed with the calculation of anharmonic and large 
amplitude fluctuation effects. 

B. Vertical Excitations 

Next we turn to the TDDFT excitation energies at the equilibrium geometry (call the 
Franck-Condon point, FCP). The TDDFT predictions with the LDA energy functional are 
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shown in Table II. The 7r-7r* transition appears at 7.5 eV excitation energy with an oscillator 
strength of / = 0.30. This energy is about 0.2 eV lower that the center of gravity of the 
observed peak. Also, comparing to the measured strength of 0.42 in this region, we infer 
that additional electronic transitions are likely to be present. Indeed, the TDDFT produces 
other states in the same region of excitation. The three additional states have a Rydberg 
character with diffuse orbitals. Only the lowest of these, which is s-like and has no angular 
nodes at large distance, is optically active. It has a transition strength of / = 0.08. Adding 
this to the it-it* strength, the total strength is only 0.04 units abelow the observed value. 
We note that our oscillator strength are quite close to those obtained by quantum chemistry 
methods, eg. ref. [5]. In principle, the LDA cannot give reliable excitation energies for the 
Rydberg states, due to an incorrect asymtotic potential. However, the present results are 
within 0.3 eV of the more refined theoretical values in the last column. It might be possible 
to improve the energies by using a GGA density functional as in ref. [16], but for the present 
study we will use the LDA and confine our attention to the it-it* transition. 

C. Coupling to the nuclear coordinates 

The nuclear coordinates that have the largest effect on the it-it* transition are the C-C 
stretch coordinate and the torsion angle coordinate. We shall consider both of these in 
the sections below. The CH 2 wagging coordinates also play a role, first in mixing the it- 
it* excitation with Rydberg excitations [17,4], and then at large torsion angles providing 
a direct path from the excited PES to the ground PES [6]. Since we ignore the Rydberg 
states and do not treat the the multidimensional character of the nuclear coordinate space 
at all, we shall neglect the wagging and other modes that couple indirectly. 

1. C-C stretch motion 

Historically, the first candidate for the source of the width of the tt-tt* transition was 
the coupling to the C-C stretch coordinate. Calculating the excited PES along the stretch 
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coordinate according to eq. (2). we find the results shown in Fig. 1. The excited PES has 
a minimum at R C c = 1-5 A with an energy 7.0 eV. The slope of the PES at the FCP is 
given by dE ex /dQ = 8 eV/A. This is all we need to estimate the width using eq. (4). Taking 
erg = 0.04 A from empirical vibrational frequency, we find a variance of oe ~ 0.3 eV for 
the absorption peak. The corresponding width is T = 0.7 eV, roughly 2/3 of the observed 
peak width. We expect that the widths from fluctuation of different coordinates add in 
quadrature. Thus, the width due to fluctuation in the C-C distance must be augmented by 
a similar width from other coordinates. 

The observed peak has an apparent asymmetry, falling off more steeply on the low-energy 
side. In principle the asymmetry can be calculated using the full expression eq. (3). We 
first solve the Schroedinger equation on the ground PES to get the wave function ^(Q), and 
then apply eq. (3) directly. The result is shown in Fig. 2. Indeed, a significant asymmetry is 
predicted, reflecting the anharmonicity in the ground PES and the curvature of the excited 
PES. However the width is still too small. We mention that the calculation of ref. [14] also 
predict an asymmetric peak. Their width includes all the vibrational modes, and comes out 
somewhat larger than the observed width. 



2. Coupling to torsional motion 

Experimentally, the strong absorption peak at 7-8 eV has a progression of strength 
at low energies, seen down to nearly 6 eV [19]. The C-C stretch PES doesn't drop to 
that low an energy, leaving the torsional mode as the most likely source of the strength. 
Combining the inertia Ig (with h 2 / Ig xs 0.0048 eV/r 2 ) of the torsional motion with empirical 
vibrational frequency ootorsion, the zero-point fluctuation variance in the torsional angle is 
given by (0 2 ) 1 / 2 = ^Jh/2I g u t orsion ~ 8°. This is rather soft, and suggests that coupling 
to this coordinate could be significant. Unfortunately, the n-n* excitation mixes with the 
Rydberg states at finite torsional angles [4,17], splitting the original ir-n* transition into 
several components. Since the Rydberg excitation energies are not accurate in the LDA, 
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we choose to ignore the splitting of the tt-ty* state, taking the diabatic state at the mean 
excitation energy associated with an excitation along the x coordinate. The diabatic excited 
PES is shown in Fig. 3. Using it in eq. (3), we obtain the strength distribution shown in 
Fig. 4. Qualitatively the measured fall-off is reproduced, but in detail the predictions are not 
as accurate as we had found for the forbidden transitions in benzene. In the region of 6.5-7.0 
eV, the theoretical strength is too high, by about a factor of three. The agreement is much 
better at lower energies, down to 6.1 eV, the lowest energy measured in ref. [19]. However, 
the theory should be more reliable the closer to the FCP, so the low energy agreement is 
probably a fortuitous cancellation of opposite-sign errors. 

A likely cause of the discrepancy just below 7 eV is the error in the TDDFT excitation 
energy. In Fig. 2, we saw that the strength around the main peak in the LDA calculation 
is somewhat lower than the observed one. This is within a typical error of the TDDFT 
for excitation energy. A overall shift of 0.2 eV in the excited state PES would improve the 
description of strength function both in the peak region and in the region as far down as 6.5 
eV. Shifting the excited state PES in that way, on the other hand, would destroy the present 
agreement below 6.5 eV. However, this large-angle region of PES is only accessible by a deep 
tunneling of the nuclear motion. It might be that the naive reflection treatment becomes 
inaccurate. For example, it might be that the tunneling in the excited state vibrational 
wave function becomes important. It may also be that thermally promoted tunneling on 
the ground PES becomes significant in the deep tunneling region, even though the thermal 
excitation probabilities are small. However, we leave it for future work to investigate these 
possibilities. 

IV. DYNAMICS ON THE EXCITED PES 

We now discuss briefly about nuclear dynamics on our excited PES. An interesting 
measurement of a dynamic time scale on the excited PES was reported in ref. [20] , made by 
a two-photon pump-probe technique. In the experiment the two photon energies were 6.2 
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eV and 4.64 eV, and the average time delay between them could be varied. The ionization 
rate was measured as a function of the delay time. The observed ionization yield peaked 
when the pulse associated with the lower energy photon was delayed by about 50 fs from the 
pulse of the more energetic photon. Clearly, the ethylene must first be excited by the 6.2 
eV photon before it can absorb the lower energy photon. The experiment was interpreted 
in terms of lifetime of the state excited by the 6.2 eV photon, with an extracted value of 
30 ± 15 fs. The ionization threshold of ethylene is 10.5 eV, so the total energy available is 
only 0.3 eV above the threshold. 

While the complete dynamics on the excited state PES are beyond the TDDFT, the 
PES we have already constructed can be used to estimate the accelerations in the short- 
time domain after excitation. According to the calculation in the previous section, the 
absorption strength at 6.2 eV comes from the torsional zero-point fluctuation the hydrogens 
rotated to a angle of about 35°. Let us assume that the ionization potential of the molecule 
is independent of the angle. Then, as long as the excitation energy is kept in the electronic 
degrees of freedom, the ionization can take place. However, the gradient in the PES will 
cause an acceleration in that coordinate, converting energy from potential to kinetic. The 
kinetic energy associated with nuclear motion cannot be used in the ionization process, 
so these accelerations give the first quenching of the ionization rate. To make a concrete 
estimate, we use the slope of the TDDFT PES to determine the acceleration. At 35°, the 
PES is still rapidly falling, with a slope given by dE ex /d6 xs 2 eV/rad 2 . Taking this as the 
force in Newton's equation, the system acquires a angular velocity given by 

d9 _ 1 dE ex 
dt ~ T e d6 

with a corresponding kinetic energy 




Evaluating these expression numerically, we find that it only takes 5 fs for the kinetic energy 
to grow to 0.5 eV. This would put the potential energy below the ionization threshold, and 
indicates that the quenching should happen very rapidly. 
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However, the fact that vibronic structure is seen in the strength function tells one that 
the motion along the torsional coordinate is not very dissipative. The system can accelerate 
to very large angle, bounce, and return to a state with no kinetic energy in the nuclear 
motion. Since we cannot say what the true irreversible rates are, the considerations here 
can only give a lower bound on the effects of a time delay. We note that the calculation of 
ref. [6] treats the large angle motion as well and obtains results consistent with the observed 
time correlation. [A 

V. CONCLUSIONS 

In this article, we aimed to elucidate how well the TDDFT works in describing electron- 
nucleus coupling. For this purpose, we analyzed the n-n* transition of ethylene molecule 
in the 6-8 eV region. The DFT and TDDFT were employed to construct the ground and 
excited PES in the LDA. A simplified treatment in the reflection approximation is made 
for the Frank-Condon factor. The TDDFT describes reasonably well the vertical excitation 
energy and the magnitude of oscillator strength. The roles of the C-C stretch motion and 
the torsional motion on the spectral line-shape were investigated. We found that the C-C 
stretch coordinate contributes substantially to the width of the transition, but the torsional 
coordinate is responsible for the low energy tail. There the transition strength is reproduced 
by about a factor of three over 4 orders of magnitude. We thus conclude that the TDDFT 
survives the first tests as a useful approximation to treat electron-nucleus couplings in excited 
states. As next steps, a more sophisticated treatment of the quantum nuclear motion would 
be required, as well as an improved treatment of the exchange-correlation potential beyond 
the LDA. 
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TABLES 



TABLE I. Ground state properties of ethylene. 

Experimental DFT Ref. [6,14] Ref. [5] 

Rcc (A) 1.339 1.335 1.347 

locc (eV) 0.201 0.19 0.21 0.22 

^torsion (eV) 0.127 0.12 0.134 0.135 



TABLE II. Calculated electronic excitations in ethylene at the FCP. The lowest 4 excitations 
are mainly configurations made by exciting a tt electron to the tt* or to a Rydberg orbital. We 



follow the convention used by ref. [13] for the x, y, z axes. 

Excited orbital to (LDA) f (LDA) Ref. [5] Ref. [21] 

s-Rydberg 6.8 0.08 7.3,0.08 6.8 

Py-Rydberg 7.3 7.5 

p z -Rydberg 7.4 7.5 

tt* 7.5 eV 0.30 8.1,0.36 7.8 

Pz-Rydberg 7.8 
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FIG. 3. Excited PES along the torsional coordinate. The solid line is the TDDFT result with 
other coordinates fixed at the FCP. The dashed line is from ref. [18]. 
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FIG. 4. Low energy absorption strength associated with the zero-point motion in the torsional 
coordinate. Solid: present theory; dashed: experiment [19]. 
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